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A dedicated effort by both the experimental and theoretical communities is crucial for achiev- 
ing a precise determination of Standard Model parameters such as the W mass (Mw)- Mw is 
measured directly at the CERN LEP2 e + e~ and the Fermilab Tevatron pp colliders, resulting in 
a precision of 5Mw/Mw = 0.03% A complementary Mw measurement is provided by the 
NuTeV collaboration which extract sin 2 #v^, and thus Mw, from the ratio of deep-inelastic 

neutral and charged-current neutrino(anti-neutrino)-Nucleon (i/N(&N)) scattering cross sections. 
However, their result differs from direct measurements performed at LEP2 and the Tevatron by 
about three standard deviations 0, S] • Possible sources for the origin of this discrepancy have been 
extensively studied in the literature (see, e. g., [6(), among them the impact of electroweak radiative 
corrections 0, H, 13 • Here we provide first (preliminary) results of a new calculation of electroweak 
0(ot) corrections with emphasis on the effects of non-zero muon and charm quark masses. We find 
non-negligible shifts in sin 2 Ow due to these mass effects but more detailed studies including de- 
tector resolution effects are needed to determine their impact on Mw as extracted by the NuTeV 
collaboration. 



I. INTRODUCTION 

The Standard Model (SM) represents the best cur- 
rent understanding of electroweak and strong interac- 
tions of elementary particles. In recent years it has 
been impressively confirmed experimentally through 
the precise determination of W and Z boson proper- 
ties at the CERN LEP and the Stanford Linear e + e~ 
colliders, and the discovery of the top quark at the 
Fermilab Tevatron pp collider. 

A precise measurement of Mw does not only pro- 
vide a further precisely known SM input parameter, 
but significantly improves the indirect limit on the 
Higgs-boson mass obtained by comparing SM predic- 
tions with electroweak precision data as illustrated in 

Fig. a 

A measurement of Mw can also be extracted from a 
measurement of the sine squared of the weak mixing 
angle, sin 2 #y^, via the well-known relation between 
the W and Z boson masses, M^ = M§(1 — sin 2 Ow)- 
The NuTeV collaboration extracts Mw from the ra- 
tio of neutral and charged-current neutrino and anti- 
neutrino cross sections @, Q. Their results differ 
from direct measurements performed at LEP2 and the 
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Tevatron by about 3cr [2|, |3[ as shown in Fig. El 

Much effort, both experimentally and theoretically, 
has gone into understanding this discrepancy. These 
efforts include studies of QCD corrections, parton dis- 
tribution functions [J, |5j, and nuclear structure (see, 
e. g., [6] for an overview). However, the impact of 
electroweak radiative corrections has not been fully 
studied yet. In the extraction of Mw from NuTeV 
data, only part of the electroweak corrections have 
been included [10]. Since then the complete calcula- 
tion of these corrections has been made available in 
the literature 0, H, 0, but a realistic, experimental 
study of their impact on the NuTeV measurement on 
Mw has not been performed yet. 

In order to remedy this situation we calculated 
the complete O(a) contribution to neutrino-nucleon 
scattering including the full muon and charm-quark 
mass dependence, which has been neglected in previ- 
ous studies. Here we present first preliminary results 
of this new calculation with emphasis on the above- 
mentioned mass effects. A detailed study, also taking 
into account more realistic detector resolution effects, 
is in progress 



II. SOME DETAILS OF THE CALCULATION 

Our calculation of the complete O(a) corrections 
to the neutral- current (NC) and charged-current (CC) 
uN(uN) scattering processes (the tree- level Feynman 
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FIG. 1: The SM prediction for Mw with dependence on 
the top-quark mass (Mt) and Higgs boson mass (Mh), 
resulting in the shaded band, is compared with the ex- 
perimental values of Mw and Mt (solid ellipse) and an 
indirect measurement from all electroweak precision data 
(dotted ellipse) Present values of Mw, and M t favor 
a relatively light SM Higgs boson, while the NuTeV value 
of M w (= 80.136 ±0.084 GeV) prefers a much higher 

Higgs boson mass. 
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FIG. 2: Direct (Tevatron, LEP2, and NuTeV) and indirect 
measurements of Mw- The NuTeV value of Mw differs 
from the world average value by about 3a [l|]. 



diagrams are shown in Fig. [3]) follows closely the treat- 
ment of 8-channel W and Z production at hadron col- 
liders of 0,0]. The O(a) corrections consist of the 
full set of electroweak one-loop diagrams and real pho- 
ton radiation from both the external charged fermion 
legs and the internal W boson in the CC process. As 
usual, they exhibit UV and IR divergences. UV diver- 
gences are canceled by including the counterterms of 
the on-shell renormalization scheme [3, EH . By ap- 
plying the two-cutoff phase-space-slicing method [16| , 
we extract the soft and collinear singularities from the 




FIG. 3: Feynman diagrams for the tree-level NC (left) and 
CC (right) isN scattering processes. 



real photonic corrections. We use fermion masses and 
a fictitious photon mass as regulators for the soft and 
collinear singularities. The photon mass dependence 
cancels in the sum of virtual and real soft-photon ra- 
diation, but mass singularities of the form log(i/m'j) 
may survive, which arise when the photon is emitted 
collinear with the charged fermion. In the case of final- 
state photon radiation, in inclusive observables these 
mass singularities cancel. However, mass singularities 
connected to initial-state photon radiation survive in 
general. These need to be absorbed in the parton 
distribution functions (PDF), which can be done in 
analogy to gluon radiation in QCD. Finally, the nu- 
merical phase space integration was done using Monte 
Carlo integration techniques based on the Vegas algo- 
rithm 

After convolution with the quark PDFs, the pre- 
dictions for the hadronic, electroweak (EW) next-to- 
leading order (NLO) cross section for vN scattering 
is obtained as follows (j = iVC, CC): 

da](E u ) = J2 j dx qi {x,Q 2 ) (da^ + da j v+s ) 

i 

+ (!•«")** 

+ J2 I dx qi (x,Q 2 ) dd-{ , (1) 

i ^ 

where the parton level cross section consists of the 
tree- level cross section, virtual, soft and collinear O(a) 
contributions (including the PDF counterterms) and 
the real hard photon radiation contribution. 



A. Fermion- mass effects 

We performed the calculation with and without in- 
cluding fermion-mass effects and are considering the 
following two cases: 

case 1: All external fermions are considered to be mass- 
less and we only keep non-zero fermion masses 
as regulators of the collinear singularities. 
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FIG. 4: Feynman diagrams for self-energy corrections to 
the NC and CC vN production processes. 



case 2: The full muon and charm-quark mass depen- 
dence is taken into account, but light external 
fermions are treated as in the first case. 

Fermion-mass effects in EW radiative corrections may 
not be numerically negligible in this process, since the 
relevant parton-level energy scale (q 2 = i) and the 
fermion masses can be of the same order of magnitude. 
This is illustrated below for the example of the W 
self-energy correction (E^) to the CC process shown 
in Fig. [H Its contribution to the one-loop corrected 
matrix element at 0(a) reads 



M 



cc 

virt 
2 



% (q 2 



M 2 W ) 



[^47 P (! - 75)^1] [^37 a ( 1 



■ 75)^2]. 
(2) 



With the renormalized W self energy E^ being de 
composed in transverse and longitudinal parts, E^ = 



[g P a - ^) ^ + Zfi-XY , one finds the following 
contribution to the NLO matrix element squared (up 

n 2 /A/f2 



to terms of 0(m 2 f /Mfr)): 
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(3) 



If we consider massless fermions for the external legs 
(case 1), the second term in Eq. ([3|) vanishes. How- 
ever, for massive fermions (case 2) the longitudinal 
two-point function contributes to the physical cross 
section. Since we work in the Feynman-' t Hooft gauge, 
we also had to include the contributions from the 
would-be Goldstone bosons, which are not explicitly 
shown here. In the s-channel W production pro- 
cess such as gauge-boson production in Drell-Yan pro- 
cesses at the Tevatron and the LHC, i is replaced with 
s, so that the second term is usually negligible. In 
t-channel deep inelastic scattering, however, fermion- 
mass effects deserve a closer investigation, especially 
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FIG. 5: Feynman diagrams for real photon radiation in 
the CC vN scattering process. Shown here are only 
exchange diagrams. 



in the small t region, which corresponds to a small 
momentum fraction x. Note that similar effects also 
arise from vertex and box corrections. 

In case of real photon radiation amplitude-level 
fermion-mass effects only arise in the CC vN scat- 
tering process (the matrix elements to real photon ra- 
diation in the NC process of cases 1 and 2 are iden- 
tical). The Feynman diagrams for the W ± exchange 
contribution are shown in Fig. and the correspond- 
ing matrix element (A4^ c ) can be written in terms of 
Z7(l)-conserved leptonic and hadronic currents as fol- 
lows (e p (k) denotes the photon polarization vector): 



M 
with 
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where jl = 1 — 75 and, 
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FIG. 7: Numerical instabilities occurr in fermion-fermion- 
photon exchange vertex corrections (green curve) owing to 
the photon propagator, while the Taylor expansion yields 
numerical stable results (red curve). 



FIG. 6: Feynman diagrams for the mixed W — <j> ex- 
change contribution to real photon radiation in the CC vN 
scattering process, and the Feynman rules for the 4 >± ff f 
coupling. 



Jm = (™1 Ml + 75)^l] [^37 P (1 - 75)^2] 



- m 4 [u 4 (l - 75)^1] [^37 P (! - 75)^2] 

- m 2 [u4Y( 1 ~ 75)^1] [^3(1 + 75)^2] 

- m 3 [U 4 7 P (l - 75)1/1] [^3(1 - 75)^2]). 



where the subscripts, I, H, HI, and IV correspond to the 
diagrams shown in Fig. [5l The fermion-mass depen- 
dence of M^ c described by vanishes when the 
contribution of the mixed — W ± exchange dia- 
grams shown in Fig. [6] is included. The only surviv- 
ing fermion-mass dependence at the amplitude-level is 
due to the ± exchange diagrams which are, however, 
suppressed by 0(m?/M^). 



B. Treatment of numerical instabilities 

It is well-known (see, e. g., Ref. [7]) that the EW 
NLO cross section to the NC process suffers from a nu- 
merical instability at small values of t owing to photon 
exchange diagrams such as shown in Fig. [4] and Fig. [3 
As a remedy of this kind of instability, we apply a Tay- 
lor expansion around small t. In Fig. [71 we illustrate 
the stability of this expansion. 

Box corrections exhibit numerical instabilities orig- 
inating from vanishing Gram determinants at small 
kinematic variables when the standard Passarino- 
Veltman reduction formalism [18[ is employed to de- 
termine the coefficient functions of vector and tensor 
four-point integrals. Especially, the crossed box con- 
tribution in the low x-region suffers from these insta- 
bilities which can be traced back to numerical unstable 
coefficients of three-point integrals. In the following 




^ »4wt ^ ^ ^ > *^ ' ^ " - 




0.0005 0.0010 0.0015 0.0020 

x (=-t/s) 




FIG. 8: The standard Passarino-Veltman reduction ex- 
hibits a numerical instability (red curve) while Eq. J5} 
provides stable results (blue curve). Shown are the three 
point functions as a function of the variable x = —t/s. 



we present a reduction formalism which yields stable 
results as illustrated in Fig. [8] for case 1 (see also, 
e. g., [21] for an alternate solution). In phase-space 
regions of small kinematic variables, the coefficients 
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of the three-point vector and tensor integral 
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can be approximated in terms of two-point functions 
as follows 
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and the scalar three-point integral reads: 
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+ 0((p 2 ) 2 ,(p 2 2 ) 2 ,( Pl -p 2 ) 2 ). (5) 

Here, 
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Pi 



B fJb1/ ...(p i ,m Q ,m i ). 



The detailed derivation of these expressions can be 
found in [l9j . In Fig. [8] we show a comparison of 
these two derivations of the Cij functions in the crit- 
ical phase space region, i. e. at small x. 



III. NUMERICAL RESULTS 

The measured value of sin 2 Ow can be extracted 
from the Paschos-Wolfenstein relation [20[ 
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(J v NC (vN vX) - (J D NC (vN vX) 
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FIG. 9: Ratio of isN scattering cross sections calculated 
with massless and massive external fermions for both the 
CC and NC production process. 



In the on- mass-shell renormalization scheme, sin 2 Ow 
is related to the masses of both W and Z bosons as 
follows: 



sin Ow = 1 



w 



M 2 



(7) 



In order to determine the impact of the higher-order 
corrections under investigation on the extracted value 
of sin 2 Ow , we use the following expression [10| : 



\^r?0 - ^- sin2 ^ + jsin 4 ^ 
A sm W - 40 • 9 n V dH NC + ° K cc) • 



40 • 2 n 
27 w 



(8) 

where the ratios of 5R v NCl SR^ C and R% are defined 
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For the numerical evaluation we use the same input 
parameters and PDF's as Ref. [7J. 

As discussed earlier, for the massless calculation 
(case 1), we neglect fermion masses whenever it is pos- 
sible, i. e. we take their masses only for regularizing 
singularities, while, in the massive case (case 2), we 
keep the muon and charm-quark masses. Note that 
following results are preliminary. As can be seen in 
Fig. [91 where we compare the vN scattering cross sec- 
tions calculated in both cases for both CC and NC 
scattering processes, fermion-mass effects are only vis- 
ible at small x and are more pronounced in the CC 
case. How this may translate into the extracted value 
of sin 2 Ow is illustrated in Tab.H More detailed stud- 
ies are under way. 



IV. CONCLUSION 

Deep-inelastic neutrino-nucleon scattering provides 
an excellent testing ground for the electroweak SM, 
complementary to e + e~ and hadronic colliders. Mea- 
surements of electroweak parameters in neutrino- 
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TABLE I: Preliminary results obtained with the 
MRST2004QED PDF set and a cut on y = -i/s > 0.12. 





Ro 


SRnc 


SR u cc 


A sin 2 0w 


Case 1 


0.30638 


0.0527 


-0.0916 


-0.0182 


Case 2 


0.31477 


0.0548 


-0.1059 


-0.0258 



nucleon scattering are not only comparable in preci- 
sion but they also probe the EW SM at many orders of 
magnitude in t = q 2 , i.e. the parton- level momentum 
transfer in these processes. 

The NuTeV collaboration used the calculation of 
Ref. [10[, which is based on a massless fermion ap- 
proximation, and did not include the entire set of 
electroweak O(a) corrections. Seventeen years later, 
a complete calculation of the O(a) corrections to 
neutrino- nucleon scattering became available [7] . In a 
follow-up paper [9j], leading higher order corrections, 
i.e. beyond one-loop, have been included as well. 
In Ref. [7[, the discussion focused on the EW input 
scheme dependence of sin 2 0w measured in neutrino- 
nucleon scattering. They concluded that the theoreti- 
cal uncertainty due to missing higher-order corrections 
has been underestimated by the NuTeV collaboration, 
and, thus, is a potential source for at least part of the 
observed discrepancy. 

In this study we focus on another potential source 
of a theoretical uncertainty, which has not been con- 



sidered before, i.e. the effects of muon and charm- 
quark masses in the calculation of electroweak correc- 
tions. We calculated the complete electroweak O(a) 
corrections to neutrino-nucleon scattering with and 
without taking into account these fermion-mass ef- 
fects. We studied their impact on sin 2 0w as extracted 
from the vN scattering cross section by the NuTeV 
collaboration. We found non-negligible differences in 
sin 2 0w when using our calculation with and without 
considering non-zero muon and charm-quark masses. 
However, a more realistic study is needed including a 
simulation of the detector resolution, for instance, to 
determine whether these effects can account for part 
of the NuTeV anomaly. Such a study is currently in 
progress. 
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